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SUMMARY 


This report describes work performed in the School 
of Mechanical Engineering at the Georgia Institute of 
Technology under NASA grant NSG-2054. The project was 
started in January of 1975 and completed in March of 
1977. Two M.S. theses [8,28], which are directly related 
to the project, have been published and it is expected 
that at least one paper will be published in a technical 
journal as a result of work performed under the grant. 

The main goal of the work has been to study theoret- 
ically the transient behavior of very low temperature heat 
pipes. A powerful technique has been developed for pre- 
dicting transient three dimensional temperature distrib- 
utions and heat fluxes for heat pipes operating under a 
wide variety of conditions. Radiating, conducting, and 
convecting external environments can be handled. In addi- 
tion, startup from the supercritical state can be accom- 
modated . 
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INTRODUCTION 




I . 

A. Statement of Problem 


Investigations into the performance of cryogenic heat 
pipes have been in progress at the Georgia Institute of 
Technology for several years. In continuation of these 
investigations the goal of this study has been to gain a 
better understanding of the transient response of cryogenic 
heat pipes and heat pipe systems. A method for predicting 
behavior of heat pipes during startup or changes in thermal 
transport has been developed. Included in the model are 
provisions for simulating startup from temperatures above 
the critical point of the working fluid. 


B. Background 

A heat pipe is a device that transfers large quanti- 
ties of heat with a relatively small temperature gradient. 
The containing envelope has an internal wick structure and 
a liquid in near equilibrium with its vapor phase. The 
phenomena of conduction and phase change are combined to 
make high thermal conductance possible. The basic theory, 
design, and operation of heat pipes are well discussed by 
Cotter [1], Chisolm [2], and others. 

As heat is added to the evaporator of a heat pipe 


,1 
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some of the fluid in the wick evaporates thus increasing the 
vapor pressure. The increased pressure drives some of the 
warm vapor to the condenser. The condensing vapor gives up 
energy which is conducted through the wall and transferred 
from the heat pipe. Due to evaporation the liquid-vapor 
interface recedes into the capillary structure and thus 
decreases the radius of curvature of the meniscus in the 
heated end. Condensation on the capillary structure at the 
cool end increases the menisus radius in the condenser. 

This difference in radii of curvature creates the driving 
force to pump the liquid back to the evaporator [!]• 

The heat pipe startup mode, described above assumes 
the vapor density at the ambient temperature is high enough 
for a continuum flow to exist and that the capillary struc- 
ture is filled with working fluid. Such startup is referred 
to as uniform startup by Cotter [3]. A second mode of start- 
up occurs if the vapor pressure is low and molecular flow 
exists at ambient temperatures. Heating produces a continuum 
in the evaporator and transition zone between the evaporator 
and condenser. Often the working fluid is in a frozen state 
and must melt before startup can occur. This mode takes 
about two and one half times longer than uniform startup 
assuming the working fluid is not frozen prior to startup. 
Non-condensable gases, when present in the vapor space, 
creates a third mode of startup. During heating the non- 
condensables are swept into the cold end of the condenser. 
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This method of startup can be accomplished very quickly or 
very slowly depending on the quantity of non-condensables 
present . 

At steady state, for some types of capillary struc- 
tures, that part of the evaporator capillary furthest from 
the condenser may be dried out. The wick geometry and 
meniscus radius determine how much of the evaporator is 
actively working. Colwell and Williams [4] suggested that 
the axial vapor Reynolds number may determine the length 
along which condensation actually occurs in the condenser. 
Assuming a saturated wick before startup the development of 
these phenomena do affect startup times. 

If the heat pipe is at a temperature above the criti- 
cal point of the working fluid before startup, all of the 
fluid will be in a single phase supercritical state. Before 
the evaporation-condensation cycle can begin the fluid must 
first be cooled below the critical temperature at the cool 
end. Capillary action will then pump the liquid along the 
wick back to the evaporator. The velocity of the fluid in 
the wick is on the order of centimeters per minute [5]. The 
entire heat pipe must be reflexed before any significant heat 
transfer can take place. 

In the design of a heat pipe system it is important to 
understand the startup of the system and the response to a 
change in power level. The parameters affecting transient 
response of a heat pipe must be determined. The goal of 
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this project was to develop a technique for predicting the 
startup and response times for various heat pipes. 

C. Literature Survey 

The performance of heat pipes in an outer space envi- 
ronment was evaluated by Kirkpatrick and Brennan [6] using 
the Applications Technology Satellite (ATS-F) . They con- 
cluded that for space applications heat pipes can be useful 
temperature control devices. Sherman and Brennan [7] per- 
formed a study of low temperature heat pipes for spacecraft 
uses. Nitrogen and oxygen were determined to be best fluids 
for operation in the cryogenic temperature range. The compos- 
ite slab wick was said to be a reliable porous media for 
present day and projected NASA applications. 

Recent work at the Georgia Institute of Technology has 
been directed to understanding the steady state performance 
of a heat pipe similar to the one in this study. Hare [8] 
performed a theoretical study of a performance of a cryogenic 
heat pipe using nitrogen as the working fluid. Polynomial 
expressions were developed for predicting the thermal prop- 
erties as a function of temperature. Mass- and heat-trans- 
port equations were combined to give capillary and sonic 
limitations on the performance of the pipe and the thermal 
resistances of each heat pipe component. A study was made 
on the effects of changing the composite slab and wick 
structure . 

Colwell [9] used a similar approach to study the effects 
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of fluid gaps in the circumferential wick. Partial dryout 
and vapor Reynolds number effect were added through a sim- 
plified approach. The effect of changing wall thickness, 
number of wick layers and mesh size was also modeled. 

A coupled heat transfer diffusion model was used by 
Saaski [26] to model a gas controlled heat pipe. The de- 
vice was 45.48 cm long with an inside radius of 1.288 cm 
and used methanol as the working fluid. Experimental re- 
sponse curves were compared with curves predicted with the 
model. Agreement is very good, both curves reaching 95 per- 
cent of steady state values in about 130 seconds. 

In a study of the dynamic behavior of heat pipes Groll, 
e_t ad^. [10] determined that small power variations around 
the operating point of a heat pipe present no special prob- 
lems. The discussion is primarily concerned with heat pipes 
in the 200 K to 500 K temperature range. A one-dimensional 
model was developed for prediction of heat pipe response and 
examples of normal startup and startup failure were presented. 
Groll concluded that the dynamic response of heat pipes can 
be adequately treated with a simple mathematical model. 

Another startup model was developed by Cotter [3] . 

Linear relationships were found for predicting uniform start- 
up and startup with non-condensable gases present. Various 
modes of heat pipe startup are described and mechanisms which 
might affect successful startup are discussed. From this 
model the response time for most high performance heat pipes 
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is between 10 and 100 seconds. 

Chato and Streckert [5] ran a series of tests with a 
water heat pipe. The device was 81,9 cm long with an inside 
diameter of 7.62 cm. Experimental values of maximum heat 
flux were in the middle of the range predicted by horizontal 

wicking tests. During normal operation, temperatures in the 

o 

vapor region did not vary by more than 1/2 C. Attempts to 
measure transients failed since the heat pipe responded 
faster than the heating and cooling system. They estimated 
the response time of the heat pipe at less than twenty seconds. 

Smirnov, et [11] developed a model for predicting 

startup of gas controlled heat pipes. Their calculations 
were in good agreement with experimental results. The dynamic 
response of buffered and unbuffered heat pipes was modeled by 
Rice and Azad [12]. Their model agreed with experimental 
data for a sodium heat pipe and for a water heat pipe. 

To operate, a heat pipe must have a source and sink. 

Any measured response must take the source and sink into 
account. Should the heat sink radiate the incoming energy 
a phenomenon referred to by Anand, e_t al_. [13] as temperature 
choking takes place. The radiating surface must adjust its 
temperature up or down to be in equilibrium with the incoming 
energy. This movement in turn forces the heat pipe to adjust 
its operating temperature. The response of such a system 
would be different from a system with a condenser bath for 
cool ing . 


I 
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Calimbas and Hulett [14] ran tests on an avionic heat 
pipe. The device used water as a working fluid. The inside 
diameter was 6.35 cm and it was 60.96 cm long. Condenser 
cooling was accomplished by a forced air heat exchanger. The 
system took fifteen minutes to adjust to a change in the 
heater level. 

A water heat pipe with a glass bead capillary structure 
was built and tested by Cosgrove, et [15]. Cooling was 

accomplished by evaporation of methanol. A steady state 
model was developed consisting of basic mass, energy and 
momentum balances. Predicted values were in close agreement 
with experimental results. After a change in power approxi- 
mately thirty-five minutes were required for the system to 
reach a steady state. 

All of the literature indicates that heat pipes have 
very small response times. Since the mass of heating and 
cooling systems is greater than that of a heat pipe they 
respond more slowly. Transient data from experiments reflect 
the response time of the heaters and cooling jackets used in 
the experiments. Any attempt to model such a system must 
include an accurate model of the heating and cooling appa- 


ratus . 
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II. THEORETICAL DEVELOPMENT 

This study is concerned with a cryogenic heat pipe 
with a composite slab wick. Geometry and physical charac- 
teristics were suggested by the sponsoring agency. Dimen- 
sions and materials are given in Table 1 and Figure 1. The 
configuration of the composite slab and wick is shown in 
Figures 2 and 3. 

The response of this heat pipe has been studied on 
analog and digital computers. A system of nodes is used to 
develop a set of time dependent equations accounting for the 
heat capacity of the system. Due to a limited capacity, the 
nodal system used for the analog model is fairly simple and 
limited to small transient. The model assumes constant 
properties and does not include any effects of fluid dynamics. 
The digital model includes a much finer nodal system and 
fluid dynamics. Property equations in the program permit the 
variation of properties with temperature. Both models pre- 
dict response and temperature profiles for the heat pipe. 

A. Description of the Digital Model 

For the digital model the heat pipe is broken into 
its three main parts: evaporator, vapor region, and condenser 

which includes the adiabatic section. Figure 4 shows the 
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Table 1, Dimensions and Materials for the Heat Pipe 
Considered in this Study 


Total length 

.9144 meters 

Length of evaporator 

.1524 meters 

Length of condenser 

.3048 meters 

Length of adiabatic section 

.4572 meters 

Outside diameter of pipe 

.635 centimeters 

Wall thickness 

.1015 centimeters 

Working fluid 

Nitrogen 

Material of pipe and capillary 
structure 

304 stainless steel 

Wick configuration 

circumferential wick 
with composite central 
slab 

Slab 

4 layers 400 mesh around 

5 layers 30 mesh screen 

Circumferential wick 

2 layers 400 mesh screen 

Fluid gap 

.003048 centimeters 


I 

I 




Circumferential Wick 



at Heat Transfer Section 
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Figure 3. Capillary Structure 
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Figure 4, Heat Pipe Model 
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parts of the heat pipe and the nodal system used for solutions. 
Basic heat transfer equations are written for each part and 
the equations are coupled at the boundaries. 

A two-dimensional model of the evaporator has been 
developed that permits temperature variation in the radial 
and circumferential directions. Symmetry is assumed such that 
only a 90° section is modeled. In this way the effects of 
partial burn-out in the circumferential wick can be studied. 

The vapor region is modeled as a lumped mass system. 

The system includes the vapor, the slab, and the innermost 
part of the wick along the entire length of the heat pipe. 

The vapor temperature is assumed to vary only with time 
having the same value in both the evaporator and condenser. 

Condenser and adiabatic sections are modeled as two 
dimensional in the radial and axial directions. This arrange- 
ment permits study of Reynolds number effects and axial 
conduction . 

Conditions at the surface of the evaporator and 
condenser can be modeled in several ways. The surface tempera- 
ture can be fixed, a known heat flux can be made to exist at 
tlie surface, or an evaporator or condenser saddle can be 
added. The condenser saddle can be cooled with a cooling 
jacket or a radiating surface. 

B. Solution Methods 

Basic heat transfer equations can be written for each 
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section of the heat pipe. The result is a system of three 
coupled differential equations: 

Evaporator : 


1 

a 





r 8r 


+ 





( 2 . 1 ) 


Vapor : 


dTy 

”vS- 3t- ' 


3T 

aF 


aT, 




C w ar 


( 2 . 2 ) 


r=r^ 


r=r- 


Condenser : 

i . 1 !IC , 

a at ^ 3z^ 


( 2 . 3 ) 


These equations are transient conduction equations and 
do not account for fluid dynamics directly. As will be shown 
in the next section, in uniform startup with a fully wetted 
wick the time for fluid velocity response is very small. In 
the case of startup from the supercritical, these effects can 
be handled by manipulation of boundary conditions in the 
nodal system. 

The equation for the vapor region is first order in 
time only. Finite difference techniques permit solving for 
the vapor temperature explicitly at each time step. The 
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equation can be written in finite difference form as: 


2TTrjj(,j,Kj^ NJE 


_ ^ V V _ 1 t w ^ ^^^n 

^n+l_^n ■ NJE Ax rj 


- T 


NKC 

I C w 


+ i L ^ fTf^ - T 

NKC Ax Cj '■^SlC-l,k V 


(2.4) 


Which can be solved explicitly 


- 2tt K., At NJE 
n+1 _ Eh ^ T 

NJE Ax m^,c^ *'^^NIE-l,j V '' 


2t: £^ K At NKC 
L 

NKC Ax in^;-Cy ^ 


- V) . (2.5) 


Note that the heat transfer into the vapor region is 
equated to the heat transfer across the inner nodes of the 
evaporator and condenser sections. 

The partial differential equations that describe heat 
transfer in the evaporator and condenser sections can be 
solved with a grid system of nodes. The curvilinear system 
is made rectangular with a change of variables. Details of 
this variable change are given in Appendix A. A heat 
balance is written for each node equating net heat gain with 
the change in temperature of that node. The resulting systems 
of equations are solved with an alternating direction implicit 
method . 


I 
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Two finite difference approximations of equation 
(2.1) are used alternately to sweep through the evaporator 
grid for each time step. These two equations are given by: 
For the first half time step 


TE^'^y^-TE^ . K 


D c r ^ ^>j = _E_ -„pn+l/2. 


+ (TE^ . , + TE^ . , - 2TE^ .) 

^y2 i,j+i 1,1-1 i,r 


K TE^ . NJE 
P 1,1 
MTE 


( E .TE*^ - 
E TEV . ^ ^ ^ ^ 

i = l 


- TC” .)C-^) 
’ Az 


(2.6) 


And for the second half time step 


T^n*^ 1 _ rpgn+1/2 

js ^ — P rTp^'^1/2 rripn+1/2 

^ -7T»r . ;nu;/2 - ^ ™i-lj - 


P c 
p p 


n+1/2 


) 


+ -2^ (TE^"^^ , + TE^"^^ . - ZTE^'^h 

^ 1.1^1 1,1-1 i,!-* 


n+1 


.n+1. 


Ay 


K TE^ . NJE „pn . f ~ i 

- .^Sh t.f i-j ■ ""‘.P'27’ 

E TE“^ . 
i=l 


( 2 . 7 ) 


Note that the equations are implicit in only one 
direction and this direction changes with each half time step. 
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The last term in each equation accounts for heat 
conducted axially along the tube. An artificial node is 
created at the evaporator end of the adiabatic section. At 
each radial location the artificial node has a temperature 
equal to the average temperature of the evaporator nodes at 
that same radial location. This node provides a constant 
temperature boundary condition for the condenser and adiabatic 
grid. A weighted fraction of the total heat conducted 
axially is subtracted from each evaporator node. 

Casting equation (2.6) in a different form gives 


AtK 1 - 1 * 1 /? ^T-i n + 1 / 7 

- c 2-) * o * 2-> ] 

2(Ax)^p c ry (Ax)y c r, 


At K „ , T At K 

( 7— E j) - ( 

2(Ax)‘^p„c„r; 


p p 1 


7-E T) (TE^ 1.1 * TE^ , ) 

2(Ay) PpCpXi^ ^ 


At K 


(1 


n 


,, — y) TE“ . 


(Ay) per. 

^p p X 


At K 


TE^ 


- (- 


2PpCp(Az) 


r) ( 


NTF 

z 

j=l 


JLlI 


-) ( 


TE^ . 


RTF 


NJE 

E TE^ 


j=l 


1.3 




(2.8) 


This equation can be used to construct a matrix of 
NIE equations for each value of j. Each equation is of the 
form : 


I 
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I 

I 

I 

I 

I 
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I 

I 

I 

I 

I 


A. + B. + C. = D. 

1 1-1, J 1 i,J 1 1 + 1, J 1 


(2.9) 


and the resulting matrix is shown below: 



This matrix is tridiagonal and can be solved with a modified 
Gaussian elimination technique [16]. This particular sub- 
routine permits any combination of boundary conditions at 
either end of the array. A matrix is constructed and solved 
for each value of j giving values at each location for the 
first half time step. The second equation (2.7) is then used 
to construct and solve a matrix of NJE equations for each 
value of i. This solution technique is known as the alter- 
nating direction implicit. It has the advantage of being 
very stable with respect to time and yet being fairly simple 
to solve numerically. 

A similar derivation in the condenser and adiabatic 
section gives : 


} 
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In the formulation 
properties are assumed constant. Expressions were developed 
by Hare . [8] which permitted evaluation of properties as a. 
function of temperature. These expressions have been modified 
to work with International units and additional expressions 
have been developed for specific heats. Appendix B lists the 
polynomials used. These expressions are used to calculate 
new properties whenever the vapor temperature has changed by 
more than a specified amount. 

Special coefficients are required for the interface 
of the pipe wall and the wick. These expressions are 
developed in Appendix C . 

The three equations were coupled at the boundaries. 
Figure 5 is a schematic diagram of the computational grid 
and boundary conditions . The temperature of the innermost 
node of both evaporator and condenser is calculated as a part 
of the vapor region and held as a constant temperature 
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of the equations the thermal 
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Figure 5. Computation Grid 
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boundary for both regions. Due to synunetry both ends of 
the evaporator grid can be assumed adiabatic. Also the 
outer surface of the adiabatic section and the lower end of 
the condenser are insulated. The vapor equation is in time 
only and therefore requires only an initial condition. 
Boundary conditions at the outer surface of the evaporator 
and condenser can be modified to better model the total 
system of which the heat pipe is a part. 

Cooling of the condenser saddle can be accomplished 
eitlier of two ways . For a radiating condenser a heat balance 
gives 

4 

^space 

which can be solved explicitly as: 


^ n*l . 

(TcSic,k-y"> 

Aj^ocAt 

^ k=NKA+l 



it 


^ ^space 

+ T ^ 






If a cooling jacket is to be used the following equation has 
to be solved: 

dl^T 
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or 
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This equation assumes a cooling jacket as shown in Figure 6. 

Both of these heat sink models are coupled to the condenser 
with a surface coefficient. 

C. Fluid Dynamic Effects 

If the startup of a heat pipe is to proceed smoothly 
fluid must be able to reach the evaporator capillary struc- 
ture at the same rate that evaporation occurs. Fluid flov/ in 
the slab can be modeled as shown in Figure 7. A one dimensional 
momentum balance can be written for the case where evaporator 
burnout has occurred: 



d 

at 





C3.1) 


If the length is assumed equal to the distance from 
the midpoint of the condenser to the midpoint of the evapora- 
tor the equation can be solved for the velocity as a function 
of time. The result is: 


2 4 
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Figure 6. Cooling Jacket Used in Model 
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V(t) 



Ky „ r L 
^5, p 


-(K 

(1-e 



) 


(3.2) 


From this expression it can be seen that the time 
required for the fluid to reach a steady state velocity in 
an initially wetted capillary structure is a function of 
viscosity, inverse permeability, and density. For the 
particular heat pipe under study the fluid velocity reaches 
better than ninety percent of its steady state value in less 
than 0.05 seconds. Although this model is no doubt over 
simplified it does show that fluid velocities do develop 
very rapidly and would present no problem to startup in this 
heat pipe. 

Velocity in the slab can be used to predict the 
capillary limited heat transfer rate with 


Q = p.VA^ (3.3) 

^cap 9 . c fg 

Heat transfer from the evaporator can be limited to this 
value by restricting the number of nodes interacting with 
the vapor region. 

In simulating startup from the supercritical an 
expression must be developed to predict the location of the 
fluid as it is pumped back towards the evaporator. Assume 
uniform and rapid condensation over the entire condenser 
surface when the temperature drops below the critical value . 
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Equation (3.1) can be solved explicitly in finite difference 
form : 


yU+l _ ^ 


2 

K - 

P£ 


(3.4) 


The length in equation (3.4) is allowed to vary with 
time as the fluid moves along the slab. The location of the 
fluid can be calculated from: 

t 

L = / V dt (3.5) 

o 

or in finite differences: 

At (3.6) 


An initial condition is needed for each equation. 

The initial velocity is assumed to be zero and the initial 
length is assumed to be one half the length of the condenser, 
measured from the midpoint of the condenser. 

This length L, the distance from the center of the 
condenser to the fluid front, is used to calculate the node 
through which the fluid is currently moving. All nodes 
containing fluid are coupled with the vapor region. Those 
nodes still without fluid are insulated from the vapor. This 
manipulation of boundary conditions, made possible through 
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the use of the tridiagonal subroutine, permits any combina- 
tion of nodes in the evaporator ' or condenser to actively 
working with or insulated from the vapor region. 

D. Logic of Solution 

The preceding equations were solved simultaneously on 
a Control Data Corporation Cyber-74 computer. A flow chart 
outlining the solution procedure is shown in Figure 8. The 
computer program is listed in Appendix D. 

After reading dimensions and initial conditions a 
subroutine is called that calculates the properties and 
coefficients for each equation. The time steps are then 
started, each step requiring one pass through the program. 
Each of the equations is solved in a subroutine. 

For each time step the new values of the vapor region 
and the heat sink temperatures are calculated first. These 
solutions are explicit and use evaporator and condenser 
values from the last time step. The new temperatures of the 
vapor and heat sink are then used as boundary conditions in 
the implicit solution of the evaporator and condenser system 
of equations . 

The same subroutine that solves for a new vapor 
temperature also calculates the position of the fluid in the 
slab and determines which nodes in the evaporator and 
condenser are currently active. Because the vapor solution 
is explicit and has little mass it tends to become unstable. 
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Figure 8, Flovv' Diagram of Computer Program 
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A discussion of stability is presented in Appendix E. This 
unstable tendency is overcome by allowing the program to 
calculate a smaller time step for the vapor region and iterate 
a sufficient number of times to complete the total time step. 

The equation for the heat sink is solved in a separate 
subroutine. This routine can handle either a radiating heat 
sink or a cooling jacket. The solution is explicit but the 
mass is sufficiently large that stability is no problem. 

Both the evaporator and condenser matrices are solved 
separately in the same subroutine. Inputs to the subroutine 
are temperatures at each node, coefficients, and boundary 
conditions. The subroutine calls the tridiagonal subroutine 
for solution of each row of the matrix. This technique is 
described in detail in references [16, 17, 18]. 

After each time step several checks are made. The 
current vapor temperature is compared to temperature used 
for calculation of property values. If the new value is 
different from the old by more than 0.1 K new coefficients 
are calculated using properties calculated at the new 
temperature. 

Another check is performed to determine if the system 
is changing rapidly. If the system is changing temperature 
at a rate less than that specified the time step is doubled 
and new coefficients are calculated. Thus several different 
time steps may be used during a single run. The response 
time of the vapor region is less than the time steps normally 
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used. This presents no problem since inherent in the solution 
technique is the assumption that the system response is 
linear across a single time step. Thus the vapor temperature 
merely tracks the temperature changes in the evaporator and 
condenser sections. However when the time step is changed 
the vapor response is solved explicitly with evaporator and 
condenser temperatures calculated for the old value of the 
time step. The vapor temperature will then lag the rest of 
the solution by one half of a time step slowing down the 
response of the entire system. To correct for this problem 
the change in the vapor temperature is always multiplied by 
the ratio of the new time step to the old time step. In 
essence this procedure amounts to predicting the temperatures 
of the evaporator and condenser nodes that interact with the 
vapor for the new time step. Since this change of time step 
only occurs when the system is changing at a very slow rate 
any error that occurs would be small. 

In startup from the supercritical the average specific 
lieat of the slab changes as the fluid moves through it. 
Whenever the fluid reaches a new node the values are calcu- 
lated again taking into account the additional fluid in the 
slab. 

Values used for thermal conductivity of the wick are 
a weighted average of the thermal conductivities of the fluid 
and the screen. Effective thermal conductivity of the screen- 
fluid combination is evaluated using the method of Williams 
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[24]. Details of this development are in Appendix F. 

Specific heats of the wick and slab are also weighted 
averages of the specific heats of the components. 

E. Description of the Analog Model 
A simplified model of the heat pipe was used for 
transient studies with an analog computer. This model 
assumed one dimensional conduction in the radial direction. 
Overall temperature drops and response times can be predicted 
with this technique for small perturbations. 

The approach used for the analog computer was to break 
the heat pipe into nodes. A heat balance was written for 
each node equating the net heat gain with temperature increase. 
The result is a system of first order differential equations 
in time with temperature the dependent variable. 

Since the computer capacity was limited only five nodes 
are used. Node 1 is composed of the outer half of the 
evaporator wall and is located at the evaporator surface. 

The inner half of the evaporator wall and the outer half of 
the wick make up node 2 which is located at the interface of 
the wick and wall. The vapor region, slab and adiabatic wall 
and wick are the components of node 3 located at the inter- 
face of wick and vapor region. Like node 2, node 5 is 
composed of one half of the wick and the inner half of the 
condenser wall and is located at the interface of condenser 
wick and wall. Node 6 is the outer half of the condenser 
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wall and like node 1 is located at the surface. Figure 9 
is a schematic diagram of the nodes and their locations. 

Equations for the six nodes are listed below. See 
Appendix G for a discussion of scaling the equations for use 
on the analog computer. 

Condenser cooling is accomplished by means of a jacket 
with a circulating fluid. An additional node, number 6 
permits the modeling of this jacket. 

The analog computer circuit for these equations is 
sliown in Figure 10. 
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Figure 9. Schematic Diagram of Analog Model 
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III. RESULTS AND COMPARISONS 

The digital model has been used to model the heat 
pipe studied by Saaski [26]. The heat pipe modeled had an 
evaporator and condenser each 15.24 cm long and a 15.0 cm 
adiabatic section. Methanol was used for the working fluid. 
The inside radius was 1.288 cm and the wall thickness was 
0.080 cm. One layer of 200 mesh screen was used as a circum- 
ferential wick. An air gap of 0.022 cm was maintained between 
the condenser surface and the heat sink for a per-unit- length 
gap conductance of approximately 0.0443 W/cm-K. The tempera- 
ture of the vapor region is non-dimensionalized by defining 

where T^^ is the initial vapor temperature and Ty^^^ is the 
steady state vapor temperature. Saaski’ s data reaches 95 
percent of steady state values in approximately 130 seconds. 
The predicted curve using the technique of this study responds 
more slowly reaching 95 percent of steady state in 160 
seconds. Both curves are shown in Figure 11. 

The digital model of the heat pipe has been used to 
study three different cryogenic heat pipe systems. Systems 
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Figure 11. Comparison of Model with Experimental Data 
from the Literature 
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are varied by changing boundary conditions on the condenser 
surface. The effect of varying heat pipe design parameters 
was studied using a model with the condenser surface tempera- 
ture fixed. Other systems modeled used a cooling jacket and 
a radiating surface as heat sinks. 

Similar studies have been performed on the constant 
condenser temperature heat pipe and the cooling jacket using 
the analog model. 

An even simpler model can be developed if the system 
is assumed to be a lumped mass. Details of this development 
for a heat pipe with a cooling jacket are given in Appendix 
H. The result is a simple exponential equation. This model 
should have a faster response than the more complex models, 
since no temperature gradient must develop. A comparison 
of digital, analog, and exponential models for startup from 
isothermal at 100 K to a heat load of 15 watts is shown in 
Figure 12. In each case the condenser heat sink is a cooling 
jacket with liquid nitrogen at 100 K as the coolant. All 
three models predict similar responses. 

The digital model of a system with a constant condenser 
surface temperature was used to study effects of varying design 
parameters. In each case the heat pipe is initially isothermal 
at 100 K. A heat flux of 15 watts is imposed on the evaporator 
surface and the heat pipe adjusts to steady state. Table 2 
lists the various cases studied. 

Figure 13 compares analog and digital response curves 
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Figure 12. Comparison of Exponential, Analog, and Digital 
Models for Heat Pipe with Cooling Jacket 
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Table 2, Parameters for Cases 

Studied 


Case 

No. 

Modification 



1 

Standard Heat Pipe as described 
Table 1 


2 

Wall thickness .2030 cm 



3 

Wall thickness .4060 cm 



4 

Circumferential wick of 4 

layers 400 

mesh 

5 

Circumferential wick of 8 

layers 400 

mesh ■ 

6 

Fluid gap .006096 cm 



7 

Fluid gap .012192 cm 



8 

Slab 2 layers 400 mesh 
2 layers 30 mesh 



9 

Slab 8 layers 400 mesh 
10 layers 30 mesh 




In each case the heat pipe is initially isothermal 
at 100 K. A heat flux of 15 watts is imposed on the 
evaporator surface. The only modification from the heat pipe 
of Table 1 is that listed. 
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for case 1. The slower response of the analog curve is due 
to the different nature of the nodal systems employed. 

The effect of changing the wall thickness is shown in 
Figure 14. Since the pipe wall is the major portion of the 
heat pipe mass, any change would have a great effect on heat 
pipe response. The pipe wall is also the largest thermal 
resistance in the heat pipe. Thus increasing the wall 
thickness requires a larger overall temperature difference. 
The time required for this temperature profile to develop 
sloivs response. 

Increasing the number of layers of circumferential 
wick also increases the steady state temperature drop from 
evaporator to condenser. This effect as well as the added 
mass slow response as shown in Figure 15. 

Figures 16 and 17 show the effect of changing the size 
of the fluid gap and the configuration of the composite slab. 
As the curves show the change in response time is very small. 

Due to the nature of the model the response takes the 
same time for any heat input. Figure 18 shows the response 
for a heat input of 5 watts and 15 watts. Also shown is the 
response of a heat pipe to a power change from 5 to 15 watts. 
Each curve is for a heat pipe with a constant temperature 
condenser. In all three cases 95 percent of total change is 
achieved at the same time. 

The thermal properties used in the model vary with 
temperature. The effect of this phenomena on response is 
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Figure 14. Effect of Changing Wall Thickness 
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Figure 15. Effect of Changing Number of Wick Layers 
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Figure 17. Effect of Changing Slab Composition 
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Figure 18. Comparison of Startup at 
Different Power Levels 
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shown in Figure 19. All three cases show an initially 
isothermal heat pipe with a constant condenser surface tempera- 
ture. A heat flux of 15 watts is imposed on the evaporator. 
Curves are shown for startup from 80, 100, and 120 K. At the 
lower temperature the specific heats are smaller and the 
response is faster. At 120 K the capillary limitation is 
less than 15 watts. The effect of capillary limitation is 
also shown in Figure 19. The development of the overall 
temperature difference is shown in Figure 20. Note that at 
120 K capillary burnout causes evaporator temperature to 
increase rapidly. 

A heat pipe system with a cooling jacket and one with 
a radiating surface are compared in Figures 21 and 22. Since 
the heat leaving the radiating surface is dependent on the 
temperature alone the response time is much greater. This 
demonstrates the phenomena of temperature choking referred to 
by Anand, e^ [13]. 

Startup from the supercritical is shown in Figure 23. 

A heat flux of 5 watts is imposed on a heat pipe initially 
isothermal at 130 K. Cooling is imposed at the same time 
with a jacket with a coolant at 100 K and a surface coeffi- 
c lent of 1000 W/m -K. The evaporator surface temperature 
climbs sharply until fluid reaches it. The temperature falls 
rapidly and the heat pipe begins working. 

Figure 24 compares wicking velocity calculated using 
equations 3.4 and 3.6 with velocities from the work of Symons 
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Figure 21. Comparison of Heat Pipe Response witli Cooling Jacket 
or Radiating Surface Heat Sink 
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I'igure 22. Response of Radiat.ing Heat Sink 






Figure 23. Startup from Supercritical 



Figure 24. Comparison of Predicted Fluid Velocity 
with Data from Literature 
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[27], For both curves the working fluid is methanol. Since 
the capillary structures are not the same the velocities do 
not agree exactly. However they are of the same order of 
magni tude . 

Wicking velocity in the composite slab is shown in 
Figure 25. The velocities are much larger than those of 
Figure 24 because of the nature of the composite slab. The 
effective inverse permeability is relatively small due to 
the presence of the 30 mesh screen yet the pumping pressure 
is relatively high due to the outer layers of 400 mesh screen. 
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Figure 25. Velocity with Position in Composite Slab 
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IV. CONCLUSIONS AND RECOMMENDATIONS 

This study has developed a technique for predicting 
transient response of heat pipes. Predicted response curves 
were compared with experimental data from the literature and 
agreement was good. Results from the digital studies have 
been compared with those from more simple analog and exponen- 
tial models and trends are found to be similar. 

The analog model provides a good method for estimating 
response for small transients. Because of the large node 
spacing and assumption of constant property values the results 
are not as accurate as the digital model and the effects of 
fluid dynamics and startup from the supercritical could not 
be included. An attempt to improve the analog model would 
have required more analog computing capability than was 
available. 

A much better technique for modeling transients is 
provided by the digital model. Thermal properties are allowed 
to vary with temperature and the effects of fluid dynamics are 
included. It was found that in accurately modeling experimen- 
tal data it was difficult to account for response times of heat- 
ing and cooling systems and the added masses of instrumentation 
and insulation. A simple approach was developed for including 
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startup from the supercritical. 

The stability of the alternating direction implicit 
solution procedure was found to be very good. (See Appendix 
E for a discussion of stability.) The program remained 
stable for time steps as large as thirty seconds. Had an 
explicit solution of the same node system been attempted, time 
steps on the order of 10'^ seconds would have been necessary 
to maintain stability. With the ADI method some oscillations 
would occur near the boundaries when the time steps were 
large and the system was changing temperature rapidly. 

In any continuation of this work many improvements 
could be made in the computer code to decrease running time. 

In addition the equations which are used to predict wick re- 
priming after dryout could be improved. 
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APPENDIX A 


TRANSFORMATION OF COORDINATE SYSTEMS 


It is convenient to transform the cylindrical 
coordinate system that describes the heat pipe to one that 
is rectangular. In cylindrical coordinates V^T is 


+ i il + 1 ^ 

r 9r r „.2 . 2 

3 r 9(J) 9 z 


Make the following substitutions 


Thus 


X = ln(rQ/rj) 


y = (p 


z = z 


3^T ^ ^ .1 
g^2 9r ^r 9x-^ 


2 

3 T 1 


1 1 9T 

“T ■ “I 9x 
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Substitution back into V T 


or 


1 3^T 1 9T 1 3T ^ 1 3^T ^ 9^T 

T dx r r r 9y 9z 


2 2 
1 9*^1 1 9 T 

T — Z "T — 2 

r 9x^ 9y^ 


9 

9"T 
9 z 


This transformation simplifies calculation and makes finite 
difference approximation of the equations easier. Figure 26 
shows the result of this transformation. 
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APPENDIX B 

THERMODYNAiMIC PROPERTY EQUATIONS 

In his thesis Hare [8] developed polynomial expressions 
for the thermal properties of nitrogen and stainless steel 
as a function of temperature. A least squares regression 
analysis program on a Wang 720 calculator was used. These 
equations have been converted to International units and are 
presented here. Expressions for the specific heats of 
nitrogen, stainless steel and aluminum were developed using 
the same technique and data from references [23,25]. 

T = temperature in degrees Kelvin 

Properties of liquid nitrogen 
Density C^g/m^) 

= -5.7784 X 10'^° T^ + 2.45356 x lO'^ T^ -3.49017 x lO'^ T^ 

+ 8.3300 X lO'"^ T^ + 2.7794 x lO'^ T^ -3.10677 x 10 T^ 
+1.31038 X 10^ T - 1.94571 x 10^ 

7 

Viscosity [N- sec/m ) 

= -2.25319 X 10'^^ T‘^ +6.541176 x lO'® T^ -5.51203 x lO'^ T^ 
+2.708099 X 10'^ T +1.0352966 x lO'^ 

Surface tension (N/m) 

= 1.026812 x 10‘^ T^ -3.919365 x 10*^ T^ + 5.631364 x lO'^ T^ 
- 3.804096 x 10'^ T + 1.1066907 x 10*^ 
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Heat of vaporization (J/Kg) 

h. = - 3.25198085 x lO'^ + 9.1831855 x lO'"^ 
fg 

- 3.4922526 x lO'^ - 1.3994675 x 10^ + 

1.9701130 X 10^ - 1.005187 x 10^ T + 2.0678544 x 10^ 

Thermal conductivity (W/m-K) 

= 3.586226 x 10‘^*^ - 1.6795344 x lO"^ + 

3.120013 X 10'^ - 2.8859648 x 10'^ + 

1.3152021 X 10'^ T - 2.182878 
Specific heat (J/Kg-K) 

Cj^ = 3.543145 X lO'"^ - 1.0894305 x 10‘^ + 

1.2632769 x 10^ - 6.500754 x 10^ T + 1.4495738 x 10^ 

Properties of 304 stainless steel 
Thermal conductivity (w/m-K) 

K = - 2.25468 x lO'"^ T^ + 9.99792 x lO"^ T + 2.2549 
P 

Specific heat (J/Kg-K) 

C = - 1.037895 X 10'^ T^ + 4.982373 T - 1.5772696 x 10^ 

P 

Specific heat of Aluminum (J/Kg-K) 

= - 8.1727423 x lO'^ T^ - 4.3008433 x lO'^ + 

8.9683338 T - 2.8647838 x 10^ 
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APPENDIX C 


WICK -WALL INTERFACE NODES 


Due to the change of variables used (see Appendix A) 
nodes are spaced exponentially through the wall and wick. 

It was found that for best results three nodes were needed 
in the wall and two in the wick with one for the interface. 
The wick nodes are smaller than the wall nodes and a scale 
factor SF was calculated so that the same value of DX could 
be used. Figure 27 shows nodes near the interface in the 
evaporator . 

To place the fourth node on the interface of the 
wick and wall an artificial radius had to be calculated. 
This artificial radius was then used for calculating the 
locations of all nodes in the wall and the scale factor for 
the nodes in the wick. 

The artificial radius is given by: 


art 


_0 

'b 


1/(NIE-1) 


(NIE-1) 

r^-Nl'ET' 


Wick nodes are smaller than wall nodes by the factor 
SF given by: 


SF = I (NIE-2)ln(rg/rj)/ln(rQ/r^^^) 


I 
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A heat balance on node 4 gives; 


ypn+ 1/2 

4 ^4‘^4 '- 2 '' ^n-^l/2_^n 


K 


CTE- 1/2 


A,, A, yrp^n+1/2 . "p A _ aAx + SFAx-, AT,y,n 

J 3 ^ i)Az (TEj^. - TE^_. * j^Az ( j 


K 


Ax + SFAx-> ,~,T.n 


TE^ . T - 2 TE^ 0 

4 ,J -1 4 ,j^ 


which can be rearranged for the subroutine to give 


K At 


-( — ^ + (1 + 
2 (Ax)^(l+SF)p^C 4 r 4 ^ 


K At 


2(Ax)‘ (l*SF)p^c^r^' 


K At 
w 


n+1/2 


2(Ax)"SF(l+SF)P4C4r^ 




-(- 


K At 

w 


-n+1/2 ^ 


,n 


2 (Ax)^SF(l+SF)p^c^r^ 


TE“ ^ y CTE'. . ., 

^ 7 r Ai.-^ r- Z ^ 4 ,J +1 


2 (Ay) 


"^4,5-1 - 2TEJ^-) ^ TEj^. 


rppH + 1 / 2 -V 
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APPENDIX D 
COMPUTER CODE 
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^ A MA I'J ( INPUT, OU rr>LfT , T Av:r5z 1 N|auT , F APE F = OU f»^UT 
1 , NAT A , r ap£(^ = ii A! A , IT S f , T APE - L I E:T ) 

'0 T .i* NN nr: F^' ( n n n . rr, n ) » t cs, •«:: no ( En , : ( ?n , f < > *w r? 3i 

i,Ar^:m,flFC(^j),AjiFr,i,Avr(t^,unAvrih,i:nAxc(6,3:), avc (b • 3c j 
1 , A ( S. ) ,i:x -: (S n : ) ,C YMf m ) *Fyr (0, C> ,CYO < f, :C ) , Axmr4 (f>nPSI2 ( 6> 

I ,OCV(l"W 

no ^NON n , rn , c n , r e. r v , i c , tp , cr c , ,n^v , ^hp ,c max , k h , nje » nic , njc 

I I L L^Fin L 10 , LNA ,cx - vr , CXC ,CYC , CV , V , T L , n , A su, ISU, XL t CC* ICRIT 

1 , A , A y , ux-, A Yc, Av, p vn sns, AC , on , c I n, Ain Aic 

1 ♦ '^'VDT , PNA 3S, UT , ,jTO , IX , <r,F f L, oxc, 0 YC ,n ,r: JA ,c RE SIS T 

1 , W T , P- , ■ F , NL » HE T A , PQ^S T Y 
PFAL NL» KL , M OOT ,>^F , NFS , kW 

w p L r : < 6 , F s ) 

6S FCPNAI (“INPu T CASE NUMf^fP**) 

PEAKS ,* F NO 
WPnE(b,7S1 

75 n E ^A r F” INPUT NT, NOT") 

PF A'M S ,* ) UT , not 
yPITc (S,M51 

as FOX NAT (••INPUT HOOT,FMflSSnTS,OrP,AX,Tl,rpfsiST**) 

PfAN (S,*I NDGT,PMASSt0IS,0IR,ap,TnCP£5lSr 
wo I r F ( 6, 1 : 3 ) 

13C F OP-^AI CTYPE 1 TO INPUT ALL OATA'M 
PEA0(5 ) INout 
r F ( INPUT .r,'E . I ) SO TO ?5C 
wpif£ isnon 

201 FOPNflT r-lN^UT NI£ ,NJ£ ,fJ IC »NJC, NJA*’) 

■RP AJ n,*) NIF ,NJF,NIC»<4JC,TJJA 
w=>ite <S,3Q:) 

39C F^=-^AT (-INPUT TSI,TEI,TvI,Trl*TRI“) 

OEA-J (S,*) T SI, TE T , T V n T Cl I T^I 
WRITE (b , urn 

40C F(nPNAT (-INPUT LW F E ,L 00 , LW FC t LNA’M 
POAj (S,»F LWP‘EiLPO,LWFC,LNA 
WPITF(6,ASa) 

450 FORNAT (-input WT,QF,PFK’L, 1ETA,P0P$TY-) 

PEAU<5,^) WT ,OP,RF,NL,OETA,pOFSTY 

c initialize values 
E T = L 
TS=TSI 

00 Z J-:,NJ£ 

00 1 1=1, NIE 

1 T E ( I , J » = T -: I 

2 CONTINUE 
TV=TVI 

no ^ m,Nic 
no 3 j=i,Njc 

3 TCatJ) = TCI 

4 CONTINUE 
r s=TRi 
XL=.15Z^ 

GO TO ISC 

?5C call OATAIN 

35C nTQ = tu 

FT=n 

wpitf, n,5C:T P'^Ass, AC ,Mc>or 

5QC [ (SX,-MACS OF M“AT SIN<^*,F11* 3,1X,-Kr, A PC* ft- , p ^ . 3 , 1 X f 

,-.pLr ^ K A Tf •• , !X , IPF 15 . b, IX ,-KG/^ EC" ) 

00 o m <un 

ATC( n = 3 , ^ 
no > j=i,nj: 

5 AT:(ii=rc(r,ji*Arp(n 

6 A T ^' ( I J A T -. ( n / FLOA nNJE > 


I 
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^ . 7 ? 

= : 

J = l 

-:s 

^ F = 1 . S ^3 c “ 1 5 
P0^jrr=,59S 
S^A ;s=i.E^?? 

L wrcO-L'^F 

HP I r ' I ^ i » NO 

40^ rc *’ 'A T ( r ? , "c AG-: r,upu^*' 
call a ^ a T d 

.iP = 9.no)7'-‘^*''^*A^*TD**44CO*lTR-TI> 

CALL OJTPUT 

c SfAPT ri-l- CVOLES 

00 IC'^J N = 1,N0T 
[rojNr=icroNT4 i 

IP (lhf:.-:o.lhf:o> go to 2 : 
call p a ^ a mt p 
OP t^- r ^ 

LWFCC -LHFC 
CO TO a: 

20 IP (A:^S ( r V -00 T\/ ) .L T . . 1 ) GO TO 30 
CALL PA>An9 

OPTV=T^ 

CO TO 

30 I F ( AOS ( 3T-S) ,L T , . [ Cl) GO TO SC 

I F (N-i 3 > -j * , u: 

40 IF ( \ P3 ( f 9 -TPO) - ,C [ 1 ?) 3 J , ac fMC 

50 iprA-^crrv-rvoi-.jcio^SfeCtAO 

55 OT=2*nr 

IF (OT-3. ) 7C, , 7. ,60 
60 OT=TC 

70 CALL PAPANTP 

OPT\/ = T\/ 

50 cy::rt40T 

00 -J. 1=1. NTE 

90 AXCONf II = (TC (I .2) -AT£ ( in *<F/ lOYC^ •2*Nje*»?*C*Are II )> 
TPO=T^ 

TVJrTV 

call VA'^09 

CALL C 0OLE9 ( AC . PC ,CC, T T , Tr. , NJC, NJA , T9 . a:^ ,C I® ,OR . QC^I 

CALi. G 91 020 (TF,AXPtAYEtCJ<£tCvT,NIE,NJ^,AXCONf 
lAr:,T7,:. TC.:.»TS.l..l..l..2.1.2t2»0.CIc.l..LWFE.LC?3) 

no :: i=i,nie 
A T 9 f I I = : . 3 
00 K J=1 iNJl 

10 A t: ( n =i r: a, j) ♦at 9 (i> 

11 AT=lI)rlT:(I)/FLOAT<NJ£) 

CALL C9IU20 < 1 C, A X C . A Y C^ CXC , C YC I C . N JC . A TC . 
lAIC,Ty/,AT'_,': T = , , 3 , 1 , 1 12 , N J A , C I C ,C JC tL WFC * L NA ) 

NJAD1 = :>IJA*1 

IF (lCJ‘>a *£0.5) GO TO 12 

CO in I'r'j 

12 call J'JTP'.JT 
ICO'JNT sS 

lOCj C0NTINJ£ 

CALL OAT A OUT 
STOP 
P NO 
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J ) 'JjTi:,'. • .'/IT V 

n I M r , ; r :) Tj t - ( ^ , ■ - ) , t c i fj , ^ r > » I ( f- , ’ : > , a y ^ j l ) , a / c < ft , ? : ) , 

1 A Y c ( n , ^ < 1 , 0 X :. ( ft , ' c ) , C Y 1 l e , T. J ^ e, ^ ^ ^ ^ ^ y ^ ^ ^ ^ J ^ ^ ^ ^ ^ ^ 

1 , a: < 01 

CO:r^J ^ i‘T • r^. n, T r, T’J.TC* T^, r-io, q:? ,Orv , ,C.I-AX aa- ,MJ£ .NIC, NJC 
: , L " : . L-^Fl , / , L J A ,nx r , u YT ♦ CXi. , r.vr , V t V , T L . iL . ASlJ , bOUtXt * CCf I CO rr 
1 , A <: . a Y : , AXC ,£^YC. A v,nv ,flS ♦ IS, :>c ,nc ,c;e ,CI C, AI F. , AIC 
1 , M , '^■lA ''S, a , flT 0 , 1 X , WF ,p , r;,0<C, D YC. CL a. JA . CFESIST 
1 . w r ,;)F , , nl . a I A , CO*.- STY 

RCAl NL.KL,^"0T,<r,MFS,KW, -iC°CL AC 
C T V = T V 

<t T V .ST . l?S) ''TV=1?S 
c viscosity ii 

1 ♦“2.7‘8:qa-o*CTv<^ia3o^'9 0S'-3 
C SU®*^Af:: T'NSION (1-/'^ 

SUr-ia= l.:^6^12':-9*C^V**4-,'Tqiq:5ej^p.^#QTy*#3^^^^3l3^i^ir-4»r^TV*»2 
1 - . 5.^ j ‘4J 9-.E - 2' CT ^ ♦ 1. 1 : ftftqr 7£ - 1 

C INVF9SS orq.r/>qTLjlY OF SLAT 
F< = o.jS 5 F'J 
C PORF 5AOIJS l'^) 

ROC'<' = :.9'_5i-0 

C OFNSIT*' OF PIPF WALL (<G/'1**2) 
r^l? 

C SOFf irti: OF P WALL (WATT S£C/KG Ki 

C=-l., 3’^9SF-2*rTV**2 + L.qp?3‘T3^CTV-l.S7 72b9 6F2 
C SP'=:ririC H£U QF SADOLF (^ATT SFC/<G <) 

CS = -d.l7?7u2''c“0*CTV**3-4.^L:^-433i-3*CTV»*2^fi.96?i33/5*CTV 

C OUTSIOc RAJTUS OF ^jpc (mj 
RO=. w^ft55 

C E^HISIVITY of H£AT SIf.'K 
E *'* = . 7 ? 

C HCAT pipe total IFNCTH (^4l 
T L = .91 4u 

C tvfiPORATn-j length r^) 

FL=.1F2^ 

C COMOFNilR length < H) 

CL = .TC 43 

C VAPpo OEnSITy C<G/m**3) 

RH(TV=3.AC 

C S^Ef.IFI.; HE-TT OF VAPOP (WATT SFC/KG Kl 
c V = 1 1 7 2 . 7 : 
o 0:=^0- ^T 

o I = tL3«U. J ♦RFynl-PE TA* ( NL-1 . C) 

PSD 1= : . : 

Ni -itjr ( iL» 

0 '^ = 

FI =FL1AT ( n 

TtPl = ALOG((PP- (FI-1, C >*A.O*PF-(FI-i,r )4H*^TA) / (R3-4.C*FI*RF-(FI-D 
C* OEi A1 1 
PSM 1= r ::RH ♦PSU'^ 

<♦3 cooriN’i: 

suha = -» ; jH 
PS'JHr. . : 

NHlzNl -1 
PO «*2I=1 , (HI 
*I =^LJAT { I ) 

Tr->irAL)G( (" o.u,.:»FI*RF. (FI-l.]1*-FTAW{a^-i4.C*FI*RF-FI*BFTAJ ) 

T|j-ir r E ♦PSHr: 

<♦2 CONTINUF 

SUH » = HS'IH 

IF (NHl .£0.3) SU*^H2 Q .0 

C 


I 



73 


c 


s .r ox , j K'*: . )Y <•; , r,rn 

^a t. .1 1 ^ M [ i— 1 
NT'- 

OTC 1 ? -MI r- ? 

N j n = N j *: - 1 

0 icn^ ‘Jic-i 


NJCMl-MJC-l 

^ ^ f 1 . V H.C/i T (fJTTM?) ) / (op ) I »• 

1 ( PL 0 T ( r“ : " ‘ ) / ( : , > F L T < 0 r : V 9 n ) 


ox 1 = ( L ( ^0 /- »nr, tjs ) ) /PL or; t <m i p*i? j 
0 Y ' zs, a ^ 3 /Fl "i:: T ( 'j j- ‘>' 1 ) 

OYC-(ALJ3(^o/-^iQr,u5i) / Float (NiCM?) 

) YC- ( tl - :l) / fl OA f ( fijcM 1 - 1 ) 

c TNEF'^fiL c OM )ir: n I TV nF ^T^r^.L-’ss stlfl (watt/m k) 

<Fz- 7 . - 

C Thtphal LJ^LlurrWITY of LIOUIO MT90GFM (WATT/f^ K) 
•^L=5.9iT?2fTF-:]*CTV*’G.;.-l.f,7 3G34LL:-7*r:T\/**L,G 

"^<^£4'>c-3*CTV»*^.:fl.3iO?Ci2ifl»C TV 

(> 3 . 1 -2 S M 

C Tw:^^A|_ C^NJJCTIVITf OF GO^-TN (WATT/f^ <) 

IL <w = <l/((OF/(2.:*->F)1*(3^;*(i<L/»<F>^OF/{?,C*pF) 

C-?..)>>?.:»<L/((PF/<2.O»\F>)*<(<t/KFM((2.0*5.F)/ 

C TH^^‘1AL C2‘nuCrrVlTY OF WICK <WATT/M <) 

n< W = < • S U A 4 - k L * S IJ '-'■'' ) / ( S U f" A f 5 U - 1 P ) 

C DFOCiry OF LIUIO NTTD03FN (KG/‘-**5) 

'^"^0^=“9.7^i4:-i,*CTV»*7^2.i4 3''^6Z*7^CTV»’^-3.L9G17C-3»CTV*'''5 
1 ♦‘^,3‘^n.F-i4*CTV*‘9f2,7 7q4F-l*CTV**3-7,i:65776l*0TV^*2 
1 ^1.51J3MF3’Crv-1.94 57lii4 
C OFNSTTY of wick (KG/‘1**3» 

^ H jW= (GJ»^A»(PC^GTY»qHOL* (l^“PD^3TY>*q)fSUMR*FM0U / (SUMA^SUHB) 

C S^FCIFir h£AT of inuiO NIT'^OG^N (WATT SFC/KG <1 

CPL=3.5*7:L3F-4*CTV»*‘+-l.^a94 33 5F-l*CTV**3>l.Z6327G9£l*CTV’*2 

1 J 7 :L r^CT V 1 . 4^95 73 S Z4 

C HCAT Or VA'^3F 17:* riPN (WATT SEC/ KGI 

HFGr-^,29iqw /45---f3»CTV^*6^9.183lAPqc-4»CTV»»S 
1 -3.4)£2 5?6 £-:*CTv*»4-1. 3 99U(S7Rcl+ 1*CTV**7 
i ►i.i/;:i3:i:«'3*crv**2-i.c.3ia7E+5»CTV 
1^2.1n7i34UE*F 
20C : F0:?>AT (4r:i5.S) 

C 5 LAO P-KOP T I -: S 

0- L3= 7 .4 7 7r«5 

OZL T = 6 . 3 1 -ir - 4 
NL A=U 
HL 1 :S 
•"A^.^Od) 

C^ = ,f>332 , 

FPCTf^wu = (‘] jr - L WFC) /rj JCMI 

•^C=»>L^O= 2 .*TL*^I»( 0 CLA»r’LA*(:A*PH 0 L*CPL*FPCTKWL + (l«-£A)*ff*C) 

1 ♦O-.L i^NL3*"r‘:^'*FM0L*CPL*F-‘CTNWL*(l .-EPI »P*C) > 

C SPECIFIC HEAT OF TM' wIC< (WATT SFC/KG <> 

CWrtd.-PJOSTvi^^^C^-POPCTvrpc'-WOLI/fd. - pQk ST Y d OfPORST V-PHOL ) 
0 W = (Cw •CUtA ♦ CP L*SUM^ T / ( A f SU'IP) 
w =» ITE ( 3, Cl : . ) r ,u PL, FF CT F'WL ,fMOL 

c niMzoGinos OF cocling jacket 
PCO=. .n7 
=» r I = . : 1 3 i 

t w = , . •: 3 1 M 

;si = d i4i<^*((CL*(‘^ro'*^-‘^ri**t) (^rc** 2 “FD*» 2 i )*c*p 

1 ♦CL* ( p : I • * 7 - - n-* r dOPL* fhol > 

IF ( 10 U.L T,: -n go to :i 

^MAC3= R iASSl 
CS = 1 . 
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11 Cn-rr.ox^r^l 


r ic - OT / _'^i b r / (K*r *n 2 . • OQ ) 

a C = o * i • c; L ♦ '* 0 / ( C - E " I S T) 

aC = ')T/ (- 

AI■=-)<■•OIS/l<^•b.^H32*:L) 
z~- :^A^:r / -0) 

A^ij = 2. 

'^S'J = £<*)^‘^iJ/“MOL 

r:^'^^K-^* "A 'j N1 - I • ( Si L z'. ^ L L A * N L £L ^ h ") l / 
; ( -'K^ X ( T L - (£L L ) / ?. n »HFG 


0 0 2C? 1 = 1, Mr' M2 

L = :u: n-i + 1 

9S1 M I) = <^nnrjs* t c>o/^por,{jsi • •('"LOfl T (L-: ) / NT:_r 2) ) ♦• 2 
00 

AY:LCi.jj=<F»0T/(nYr:**2^^*c*?.*9Si?(in 

CY“(I,JI=AXl(I,]> 

ICC 0 Y M I , J ) = A Y £ ( I , J > 

20C ooNfiMu: 

SF::.5*Nr:'"'2*AL0G<S05T (‘^Sl2(Nir'^2> I / ^I»/flLCG(«-C/ P^OGUS) 

p r?='^r^*2 

PNI il? = ( - OGUS* ( = 0/f=^'iOG'JS) - M2/2 J ) • *2 

^si?{Ni-:Mi) = (ci*(5CPyrcsi2(Ni£i2n/on»*,=ii**2 
X n^H - ( ( ?5 I ? < M - M? ) -9 u;*^0 2' ♦ PHO>**cw + (^M ^1 2- PS12 < Ml £M2) ) ♦R*C) / 

1 (9Mrnc -PK‘<"'o?> 

AV = Tf^>J*U*2.*L3Y£^c“t/(nX'i*<PMOV*TL*(3.141f»PI2- 
1 ?.*^I^SlPTHX ) *CVt ’K=SL A0>T L* (Fimn2-PI2)* 

2 3 .141S^C’^) *SF ^ 

.5V=TK’^^ilT«2.*<s,2q7»OYC/fCXC*('^HOV*TL*(3.141G*PI2- 
12 . • ILoTHO • c;f MCP$L At- f TL* (FINMn2-P 12 ) • 

2 3 * 141 CW) • St^ ) 

IX^lNTlAV^NJLfJV^NjO^lC 
A V =AV/FL0A T t IX I 
nV:^PV/FLCAT { IX > 

GV = TKW«jY-*:iL/<OxE*Sn 
‘^PIT£(?,2jC0 1 0’MSS,CW,^CPSLA^tCS 

IOC : J FOP1AT 

no J-l,NJ? 

.'\x£C:<I£12,J)=K'r*nT/(nxF**2*xnSH*(l.4‘SFi*=’GI2(MlEM2l> 

CX£( ^4 :£>!?, J)^W^«nT/(PXr»*’*xrGH*(l,*S=^)^9SI2(NlFr-2) *SF > 
AY:(•^4TP>^?,J)=KF•uT/(i)X£••?•X^SH»2.*PGT?tM=:^'2n 

r.Yc , J) = Av-(nIph2,j) 

C Xt (NI El l,J» = rK4*OT/(DXE**?»SF**2*C'r4»CHOW*2»RSr2(NlEMin 
AX£(Mi:n*J)=CXE(NlE‘-l»J1 

AY:(NlE'n,J» = T<W*0T/<jYr ? • P^OW* C W* 2 . ♦ P I 2) 

CYE<n:£'^1,J)=Av-(NI"mi,j) 

AXE(n:-,J)=‘*>ST?<MI£M1)/i. I2*CXMME^1,J) 
gx£ (n: £, J ) = • xE (v: : , j> 

AY EtME,J» = TKW*0!/(OY PM OW»CW*2.»PI2) 

30C CY«= (M £, J> rAVf (Nlf ,J> 

no 50.: l = :,Mcy 2 

no j = i,^fj: 

Axr. (:,j)=KFM)T/<oxc*‘2»p*c*?.*Rsi?(in 

ay: (I, j)=‘<F*iM/invc*»2*^*r*2,> 

GXO ( I , J) = Axr n , J) 

4C c c Yc a , j» = ayl ( I , J) 


I 
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I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 



10 •'■).: Jr i^rijr 

ft : u i ■; " 1 , Ji = T < IT / < :jr^ • • n>j*2 , ) 

ft ( N! C 1? , J) =K J 1 / ( n Y C* XUSH*2 . ) 

n ^ c rc i; u . j I = ft vt: ( n i o i ? , j ) 
c Yonn o*n t j) - ft yi:< Nic*n . j i 

cxc(Ni:i?»jj=T<.v*nT/(rxc**?*y^SH’<:.*'^'^)^«sie(NicH?MSF) 
cxcnidi, j) = T<*j*[iT/(nx(:**?*^p«*?*ow‘rMnw*2*fr,i2(vic^<i» > 
ft xo ( ii cn , j) = cxc t N I CMi ,j ) 

CYO(NIOf J> = TKW*nT/ (nvc*» ?*omow^CW*2 , \ 

ft Y C (*l I 0, J > >C YL ('Hr , J ) 

ftX0<‘4lC,J)=PSI2tNT£Hl)/i^I2*CXC(NIC^l^J> 

60C CXC(^:.:, j) rftvrrjic ,J) 

?r ru^ ) 

^ NO 
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I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 



"'.’I i ->0’i r : N V ^ 

^ M ^' . I 'M r - ( f , ’ : ► , r<‘ I ^ , A ^ , A V ( b ,3 c ) f 6x c ( ^ 1 1 j ) t 
1 avo (-3 . , cx M e, i : > ,rv.' (6 • ^ : 1 ,c>cc ( 6 1 ,r Yc (e .30 ) • axcoNi-b > 

1 .3:v(i-i 

coMiT^i ' j iTs,r. r',T",Tv,Tr,T-,CT^tQ==iO'V.r:f:R,ciMx,Mi',MjE*NrG.Njc 
l.L-^P~,LHr.%L^^O,LNA ,nx ^,rYF,ryr,rYr. ,GV, V, TL.FL , ASU,RSU,XL,CC. ICRIT 
1 , AX - , Av- , .^YG , A V,;, A V . G V . A S . “>■ T- , £ C , G C , C I C t C I C , A 1 1 , A I C 

1 . OT , 3'i A , w , m , IX , xp , c'. r , r;<r , n y r , tL , ruA , C^^GSIS T 

1 f wr, J^’ ,-'r ,.JL . -"TA ,Pn-STY 
PfAL *IL, <L , '^nc r, <F , hFG ♦< W 
GO TO ^ J : 

X L n = G J' • ( 1/ f 1 MtiS (Or^: V) / G^'AX> ) ♦►I 
XLNA=(MiO-NJA)*OF\//0‘^AX + NjAfl 
L ^1= r ^ T ( X L ^ PJT ( ?. • ( XL TMT (XL 00 ) ) ) 

TF (L jO, r. T , ( M J 1-- f 1 ) ) LF0=NJr4-l 

LN A = I JT ( XL NA ) TUT ( 2. • ( XL NA-IMT ( XLN A) ) ) 

rr ( lna ,le . ( ujA ►!) ) GO to m 

IF (LO A , GT . ( NJCf I ) ) GC TO 'fTi 
GO TO ?:g 
30 LNA^NJO^l 

G 0 T 0 e : C 
IOC UN'ArrUAf2 

?0C NJAo:=NJA*l 

su nc = o. G 

00 3tt J- JAPl f NJC 
SU ^TO-SUY rc^ TO fUlCYl , J ) 

SOOTC=SiJ''TO/FL OAT (NJC-MJA) 

IF (GU.-irC.GT, 12 G) GO TO 5 G 
IF (ir,^^. I r .LT . 1) GO TO 5G 

0 0 ^0 1^1 , i: 

V=nT/iJK.»(ASU-0SU*V-’XL-\/*»2)/XL>V 
xL=KL>v*QT/::: c. 

29 CONTIOiJ-: 

IF (XL. GT . (TL-CL/2 . ») ICkIT = C 

LOF=IuT(FLOAT( NJC -NJA) / 2)>IVT (<L/jVO)>IrJT (2*(yL/CYC-INT(XL/0YC>») 

1 

LWFC=NJC-L0F 
IF (lWFC-NJ A) 3i) .73 • 3C 

30 LWFf; = NJA 
GO TO 4G 

3S IF(lwFC-I) u:,t.G.U5 

4Q U'^FC=1 

XL0)=(XL-TL*nL/3+EL»/tL*NJ? 

LO0-rJT(XLOQ)*lNT<2'^(XL'>C-rNT<XLn0)))+l 

1 F (L°0 ,G T . (N Jl f L) ) LOO=N f 1 

4S GO TO 9G 

50 Lk^Fr.r'slJG 

L'’0=1 
ir^i r = i 

55 continue 

IF (L > GO TO iCul 

Q-'/r,. 

OCVd) =.5*CV 
or‘»/(UJE) -.5TV 
Nj-n = NJE-i 

00 uui 

99 H OFV IJ) =0 V 

1^0-: 

on 1 : . : Jr 1 , vj' 

IF (no.cj.: ) oo ro ic«c 

0- 7 = n^\/( Jl • ( Tr CNTE-t , J) - Tr ^^.IT = , 

IF (oev.lt.o^axi go to i::o 

1 “>o-i 


■ 


I 
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rc^iT^i 
L'^3- J- 1 

lOCC nrn:NM- 
IG G l n n:iT i'4n\ 

A V *: I ■ - 3 . L 

A V ' T "■ V = 7 . 

A \/ 1 T G - J . 

A VE ) , c 
L = L n-1 

t wpcr>: - lviF“> 1 

E AC^ I 7 - - FlOu T ( L'^*0 M 1-L NT - Di f 1 ) 
nflCTWE=*^LOAT{LNAMl-LWFr,Pt*l» 

(U'" UL £ . LWFE P; ) r,n TO lu 
no i: j=!_’/<ff i->i , Lsn^'i 

fl VE rev =A \/ - TEV* T£ (N IP , J ) 

1? A V£ T: ^ AVF I £ ♦ T : (MI E -1 , J) 

AV£T' 7=.1V - TEV/FACT WE 
AV~T-=a;F. TE/‘ACTIV£ 

If (MJF-L'^O^n) 111 11. 12 

11 f ArjivE^tACTWr-.f 

12 IF (’“LVFEPDlSil W13 

IJ £ACTI;E-£ACT WF-,5 

Gn TO IS 

lU FACriV£-WV 

IS I F (LM^ .L E .LWFCPl) GO TO 2U 

IS 00 e: 

A ^ETC^ -A V - rC‘7^ TO ( ^’ IC I J ) 

2C AV"FCz A'/f TC^ rCTMTC -1 t J) 

AV'ET07=AV'. fCV/CACTIVE 
AVETO" Wr rC/CACTIVE 
IF<mj:-l^:a«1 )Zl,?l,22 

21 CACnvfrrACT IVE-, 5 

22 I F C?-LWFC^n 2S ,23 . 23 

23 CACTI7E^CACTIVF..S 
GO TO ?S 

24 0ACTIV£=C.: 

25 rvo=rv 

00 2S T^l fix 

2S TVrAVEAOTIVE^IAVETF-AVLTrv-TV+’TVOl^TV 
1 ♦ OV ♦OAGT IVF. • ( A\/J^ TC- T VO- !';♦ TVO ) 

TV=(TV-TV0»'01/OTOi>TVC 
0 TO-lU 

W9IT£(3,4^>L^r,LMA,UWFCfV,XL 
49 FC^'iAT r «I'i,2Eie.5» 

807 FO^ 1AT ( *F 15, 8) 

3 TliPN 
ENQ 



t'l r* 


, r T , T »;uc . NJA , E 


SJ'^T'JTIN' COOL:'^ 

>nN T (P, ■»:) 
r = : . ; 

fijft n = MJr, ♦! 

no :: j = i , mjc 

SUM T =*J'M r T < : » J) - T t? 

on- = a« sj^ T / ^^LC(w ( Njc-‘Jja ) 

TP = 1 * ( ♦!=> 

5- T‘J^N 
£N >1 



JIJTI M-- "> J I n?.) ( T , C X. /\ Y ,r V ,r;v^ fi T ,N j, I AL 1 AI 1 , IM , 1 , fi JM , 

1 Q U 'Ml , L n •!. L J1 tL ,r,J A ,r: II , r: Jl , LV^c* ,L = 01 

D I '■Mh . roM T ( s , r » , ri( M ) , v( ( 7 M M M ) .r ( ^ : ) , aji <s) i ax ( ^ 

1 A Y (q , 3 ,rx\r^,M) , CY ( fj ) , A( M N C ( 3M , AX IAl< b) 

Ml '11 =11-1 
M J‘11 = ^r J-1 

c HA<r r sx— 

c piC‘<' nij'iTA ?y coNTiriOM 

ir (lji-2) 

1^* 00 I*” : = 3Mii*n 

IS 0 ( n = A ji 1 1 1 

00 TO 1 } . '; 

ie> 00 17 I=2,NI^l 

1 ? n(I> = (U-2.*AY{T,in*T(I, !»♦?•♦ AY ^I,U’'Trif?)^T(I*lMAKlAL(n 

TM TO t:c 

133 iF(i-‘ijAi:::,20M3T’c 

?0C 

A l;>i=i: . 3 
GO ID 4 : c 

30C LT‘J = lI1 
A IN=AI 1 

4C0 oontin'j: 

IF(I-LW") S]J ,8p: , soc 

soo I (1-L “^o> , 6 JO 

60C LI0=2 

AI0=C . 0 
GO TO 0 ) D 
70C LI3 = LM 

a:o=ai^ 

800 oor^riNU’ 

00 1? j3 I = 1MJT 
A (I)=AX(IM) 

10 03 C(I>=CX(T,U 

call - 0AC-1C ( A ,G* OiNl , E ♦ F , W ,L [N,L I 0| 

1AIN,AIO,0[1.0IH) 

21 no 4'j J=?,NJ01 

OC 23 T=2,MIMi 

2s rMn = :Y(i,j)*r(i,j*i)fAY(r,ji*T(i,j-i) 

1 ♦(l.-AY(I,J)-C^(IiJ))'‘T(T,J)4^T{T*jT’AXIAL(n 



00 ■?: i-i 

,NT 

30 

M I ,J- l) = 

Ail) 

IQl 

IF ( J-N JA ) 2C1 »2C-1, 3C1 

20 1 

LTN=2 
A IN=: . c 
GO TO 43 1 
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L TN = LI 1 
A T N = A I t 
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CO U INU- 



tF (J-cWF ) 

b 3 1 , f - U 5 0 1 

531 

IF(j-l»'^0» 

^Cl.tMltbOl 

631 

L 10 = 2 
A 10 = J . } 

GO TO 1.1 


701 

L TO^LIH 
A TO = AI A 


801 

C C N T [ OJ‘ 



no 1^.1 I 

= 1 ♦ N I 


fl ( T » = A X f T 

. J» 

IOC* 

c ( n = 'X ( I 

,J> 


OALL (<*,'•., n, ‘IT F ,L IN, LiOt 


IAIN, A [0, on. OI'I) 
GONT IMU: 
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■)c S'* 

U4 D(n = (:,-?.*Gva,Nj))*nr,‘jj) + >.*AYn,Nj)M(i,Nj-i) 

1 ♦ r ( I , *4 j ) • Y T L ( n 

ur 00 50 1^1, NX 

31 T ( i.;u-: )='A{i) 

ir inj-la' » 6c ^. » 5:? 

53 2 I <’U-L <0) 7 ■ 2 , 2 

63? LI0=2 

A ro=: . 0 
GO TO i]2 
7C2 l:o=lh 

AI0=AIH 
fl02 CONTINU- 
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10C2 CO^CXII.NJ) 

call ^.riAC0£ (A,C,n,MtE:,F,W,LIl,LIO. 
iAii,fin» oil . oiM) 

nn 4*^ T-’itNI 

UH T(I,r^JJ-W(I) 

c ha<l j sw:^p 

IP (LI1-’ ) 5 1, 51 , S3 

49 0 0 5 :> J - 1 < N J 

50 n(j>=Ai: 

GO TO 50 

51 no 52 

5? 0(J)“(l.-2.*AY(l,J)l*T(l,J)f2.«flYtl.J>*Tr?,J)-CIl*AIl 

GO TO 50 

53 00 54 J=?,NJM1 

132 TP(J-OJ\) 202,20?, 3C2 

202 LI0=2 

Cl. 9 = 0.: 

GO TO 432 
30? LI9=LII 

CTN=CI1 

40 2 CC9II9(/£ 

0(J»-(l.-:.^Ay(l,J)-2.^CIN»*T(l,j)»2.*AX(l,J)*TC2*J> 
!♦ (2.*CIN)*0I1*T(i,J)*gxIAL(1) 

54 IF (LIN. £G. 1) 0 f JJ =AI1 

nn ic: 3 j=unj 

A ( J)=AV( 1 , J) 

ICC 3 c (JJ=0Y(1 , J) 

55 call >»0AC-*t <A,C*n,9J.L,P,W,LJlfLJ^f 

lA Jill) J1 ,njH) 

56 no 7c 1 = 2 , NIm: 
no Gv J=2,0J“1 
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1 +ll.-V<(I,J)-CX(I,J)lM(I,j>+T(I,J)*AXlAUm 

no ss )=i,sj 

65 T(I-l,J)=WlJ) 

no i:: 4 j=i,nj 

A ( J) = A I . J) 

1CC4 c < J)=CY( I . J) 

call r OA C H£ ( a , C * 0 , 9 j , - , F , W , L j I * L j H , 

1 A J1 (I) , IJH.CJl ,0JM) 

70 CniJTirjij: 

C PTC< 0 OTTO 4 ‘n!jNOA*.?Y r. ONOlriON 
no 75 J=ltN'J 
IF ( ) "3, 7’ ,-^2 

72 IF( J-cl31 /4,7,^,73 

73 n(Jl = (l.-'.*Gx(NT,J))*T(^I,J)+2*AX(MI,J)M(NI-l,J) 

1 ♦Tirjr,J)*AXlAL(Nn 

GO TO 75 
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I 

■ TIJ)=At^ 

* 7S CO‘JTt>»U-: 

no 

6') T rrJI -1 , J) - W( j) 

■ no j=i,Nj 

V 

1305 c ( J» =nv ^Nl ♦ j> 

I CALL (a ,c. o*MJ ,E .F,W,LJl^LJ^if 

1 A Jl( Nl) , A JM, OJ 1, OJM) 
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I 
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S U ^ ^O'l T : fi , y ^ A ' < A , c , , - 1 , - , - , w , L 1 , L , n , Ah' , 

nt ♦ VM 

.n "»■ N'l r IN A < ,c ( Vl)) , 0( mp, , c ( ^P) , c ( mD) , W {'M'.) 

di-:} 1 ,2 

1 iil)') 

r u)=Ai 

GO ro ^ 

? ^ ( n = i , 

< n Ai 
GO ro 

3 r 1 1 ) - A 1 / ( A : - 1 , ) 

^m=01/(i.-AD 

on T 1 - 7 , r-* 

^ ( .1) = 0 ( 1 ) / J F N 

IP as-no,7,8 

6 H ( 'lOJ = A ^ 

GO TO 0 

7 W ( o ) = ( ( '1 1 ) ♦ A M ) / < 1 , - r ( M M ) J 

GO TO ^ 

8 >1 C 0 ) - < F < M *i ) A h! ) / ( ( 1 . f A ^, ) / a M -E ( M M ) ) 

q no i: 

Ij W ( M) r- M ) ♦W ( 1 ) f F (M) 

tu=>:j 

“NO 
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r-o J f ! our^'jT 

n ) t rr { Ps, i: ) , ^ , A Y' (s ) . Axc<G, 10 » 

1 • flxr.(o , ^0 ) , CX- 3" ) , CY- (f' , V > t Axcn^| (o> 

O"*i0'^ *. T . T"; ,'. T^, Tf , r \/* T r , T C?, Iv, ij = ,r,r\/ , OCO tCMax , N lii ,MJE * NIC , tIJC 
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00 ^ I=:,NIC 
N=NIC-I*1 

Z uoijr (1,7Q| ( TO (N,J> t J=1 *NJC) 
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S r n - -1 ft T n 
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1 “0NAX“,F8.1,3X,"htATTS**) 
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1 f LWFT . L'^fCC ,L'n ,lNA fCV f ,C YC » CXC ♦r.YC , CV , V • TL • £L . ASU, HSUtXL ♦ CC, IGRIT 
l,AX-,/iy-,AXCiAYC,AV,nVfAS,P^,AC,nCiCIF,UCfAirfAIC 
1 . T , A ■'iC, ,)T ,() f 0,iy , <c , F , c, OXC, 0 YC, CL , fJJA ,C RFSI5 r 
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A:M > , 7b ) L4F: . LP ^ t LWFC , LNA , ir.RI T 

P - AM ( ^ . r.4 T , , °F ,NL ^ li T A , PO^STY 

p"A )(?,;: ^» rs, GS£ 

")0 1 : 1= 1 , NIC 

1? AT ( e . tTc n , j) , j = i ♦Njc) 

5£A,)(^,7r • ) rv,V.XL 

10 ^: r=iiNic 

20 »*A0 (a,t4,3) I TCI I ,J) * J = ltNJCI 

PE A')( ’owO T^OCRiQP 
59 Fr>-'‘>MAT m ■! ) 

75 FCn-^AT 

b: F0-MAr<3-:i^^5^pe,,2,2ElC.5) 

10 0 FO^TAT<?F‘^.3) 
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4?C rc- -rar 3} 

5 jC FC^OAT MF 3.3 ) 

=>r Tu**^o 
f N J 


I 


I 


t 


86 


APPENDIX E 

DISCUSSION OF STABILITY AND ACCURACY 
OF NUMERICAL TECHNIQUES 

The technique used for solution of the vapor tempera- 
ture is essentially a one dimensional heat balance of the 
form : 


■ A. 

1 1 


+ B. 

1 + 1 1 


pH 

‘i-1 


(l-A.+B. ) 
^ 11 ^ 


T. 

1 


n 


For stability it is necessary that [16,17,18] 


(A. + B.) < 1 
^ 1 1 -^ — 


All elements of A and B are fixed except the time step. 

Thus the time step must be adjusted to maintain stability. 

In solving for the vapor temperature the coefficients 
A and B are calculated. New values are then defined by 


A. 


A. = 


i ■ (a~B.T 
1 


Bi 

^i ■ (A'.+B.) 

1 


i 


I 
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The solution is then iterated (A + B) times to complete 
one time step. 

The alternating direction implicit method used in 

2 2 2 

this study has an accuracy of 0(At ,Ax ,Ay ) [16]. 

According to a von Neumann stability analysis the 
method is unconditionally stable [16,18]. However in 
practice this method, although it does permit larger time 
steps than standard explicit method, will become unstable 
if time steps used are too large. 
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APPENDIX F 

CALCULATION OF THERMAL CONDUCTIVITY FOR THE WICK 

Since the wick is a composite of screen and fluid 
gaps, a method must be developed for estimating its thermal 
conductivity , 

The thermal conductivity of the fluid gaps is given 
by the polynomial expression given in Appendix B. 

The thermal conductivity of the screen was calculated 
using the method of Williams [ 24 ]. 


K = K, 




[2 


1 

IT 

P 


^f 

2r, 


- 2 ] 




(^) ( * 

P 


3^777 


1 ) 


+ 




■3-TF7 



The overall thermal conductivity is then a weighted 
average of the values for the screen and fluid. Weighting 
factors are given by 


I 
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I screen = Z In [ (r„- (4r .+6) (N- 1) ) / (r -4Nr^- (N- 1) 3) ] 
N = 1 

NL - 1 

"o fluid = E In [ (rp-4Nrc- (N-l)B)/(rg-(4r£+6)N)] 

N=1 

Total thermal conductivity is calculated to be 

^ fluid) 
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I 

I 

I 

\ 

I 

I 


I 
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APPENDIX G 
ANALOG SCALING 

Due to the nature o£ the analog computer equations 
must be scaled before they can be solved. This scaling is 
accomplished by making a change of variable. Each variable 
in the unsealed equation is non-dimensionalized by some 
maximum value. For example the unsealed form of equation 
4.1 is: 


dT 


1 


4itLt- K 


dt p c v p In (r„/ rp) ^"^2 "^1^ ^ p c v „ 
p p pE ^ 0 B P P pE 


2Q, 


Make the following changes of variable: 


T-T . 
min 

T -'T . 
max min 


Qe 


* 


Qe 

Q — 

^max 


t* 


t 

t 

max 


gives 



(T -T . ) dT * 
^ max min^ 1 

t 

max 


*E '"p 

PpVpE 



-Tj*J 


+ 


20^* Q 

^max 

p c V c 
p p pE 


which can be simplified to: 


dT * 4tt 2.- K t 

1 _ E p max ^rj, * T’ *-v 

In (rQTrgT ^^2 ^ 


dT^ 


p c V c 
P P pE 


2 Q t 

'^max max 

p c V _ [ T -T i 

p p pE max min 


Qe* 


Similar scaling for each equation gives a new system of 
equations that differs from the first in only the time 
scaling term. 
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with c 


Define 


to give 


and 


APPENDIX H 

DERIVATION OF EXPONENTIAL MODEL 

Assume entire system approximated as lumped mass 

Then 

nipCf (Tj-T) 


jmposite specific heat 


"'t^t 3T ^IN 


T* - (T-T,). Equation becomes 


dT* 

”h-^T dt ^IN ■ mfCfT* 


This first order differential equation can be solved 


"'f^f 

Qjn 'V^t 

T* = (1-e ) 

mfCf 


Qout m^c^T* 
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